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THE cx)(x)-EQUATION IN RIEMANNIAN VECTOR 

FIELDS 


THOMAS BIESKE 


Abstract. We employ Riemannian jets which are adapted to the 
Riemannian geometry to obtain the existence-uniqueness of viscos¬ 
ity solutions to the oo(a;)-Laplace equation in Riemannian vector 
fields. Due to the differences between Euclidean jets and Riemann¬ 
ian jets, the Euclidean method of proof is not valid in this envi¬ 
ronment. 


1. Introduction 


Recently, the p(x)-Laplace equation and its limit equation, the cxd(x)- 
Laplace equation, have been the focus of much attention as a tool for 
exploring applications such as image restoration jl] and electrorheolog- 
ical fluid flow [7]. Linqvist and Luukari [6] recently proved existence- 
uniqueness of viscosity solutions to the cxD(x)-Laplace equation in (Eu¬ 
clidean) M"’. However, this proof is not valid in general Carnot- 
Caratheodory spaces, such as Riemannian vector fields, because it re¬ 
lies on two important Euclidean properties, namely that the so-called 
viscosity penalty function is the square of the intrinsic distance and 
that the two first-order jet elements derived from the penalty function 
are equal. (These two phenomena are discussed more below.) The 
main result of this paper is that the lack of these phenomena in Rie¬ 
mannian vector fields can be overcome to produce existence-uniqueness 
of viscosity solutions in this environment. In particular, we prove the 
following theorem: 

Main Theorem. Let Q be a bounded domain in ML with Riemannian 
vector fields and let f : dfl M be a (Riemannian) Lipschitz function. 
Then the Dirichlet problem 
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has a unique viscosity solution u. 

In Section 2, we review the main results and definitions from Rie- 
mannian vector fields. Section 3 is dedicated to existence-uniqueness 
of viscosity solutions and Section 4 details further properties of the 
viscosity solutions. 

2. Riemannian Vector Fields 

2.1. The Environment. To create a Riemannian space, we begin 
with and replace the Euclidean vector fields , dx„} with 

an arbitrary collection of orthonormal vector fields or frame 

X = {Vi,V 2 ,...,Vj 

consisting of n linearly independent smooth vector fields with the re¬ 
lation 

^ d 

Xi{x) = 

j=l 1 

for some choice of smooth functions aij{x). Denote by A(a;) the matrix 
whose (hj)-entry is aij{x). We always assume that det(A(a:)) 7 ^ 0 in 
M”. 

The distance between points x and y, denoted d{x,y), is defined as 
the infimum of lengths of curves that join x and y with the additional 
requirement that the curves’ tangent vectors lie in the span of the 
Xfs. Using this distance, M” with this frame is a metric space and, 
unlike an arbitrary Carnot-Caratheodory space, this distance is locally 
comparable to Euclidean distance. We will discuss the importance of 
this fact below. 

The natural gradient is the vector 

DxU={X,{u),X 2 {u),...,Xn{u)) 

and the natural second derivative is the nxn not necessarily symmetric 
matrix with entries Xi{Xj{u)). Because of the lack of symmetry, we 
introduce the symmetrized second-order derivative matrix with respect 
to this frame, given by 

{Dlu)' = \(Dlu + (Dlu)‘). 

We can define function spaces and the Sobolev spaces etc 

with respect to this frame in the usual way. 

We may also define the cxD-Laplace operator 

Ax,ooR = {{Dluf Dxu, Dxu). 
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This operator is the “limit” operator of the p-Laplace operator (for 
2 < p < oo), which is given by 

= \\Dxu\\^~^Axu + (p - 2)\\Dxu\\^~^Ax,ooU 

= divx {\\Dxu\\^~‘^Dxu) 

where the divergence is taken with respect to the frame X. 

Following [6], we generalize these operators by replacing the constant 
p with an appropriate fnnction p(x) G (7^ fl W^’°° and scalar fc > 1 to 
obtain the p(x)-Laplace operator 

Ax,p(.)W = + 

+ \\Dxuf^^^^-\Dxu, Dxkp{x)) In \\Dxu\\ 

= div^ {\\Dxu\\^^^^^-^Dxu). 


The corresponding eqnation Ax,p{x)U = 0 is the Enler-Lagrange eqna- 
tion associated to the energy fnnctional 


\\Dxu\\^P^^^ 

kp{x) 


dx 


Allowing k —)■ oo, one has the tool for analysis of the extremal problem 

min max 11 Dxu 11 

U X 

Letting k ^ oo, we have Ax^p(x)U ‘^x,oo(x)U where 

Ax,oc{x)U = Ax,ocU+ \\Dxuf {Dxu, Dxlnp{x)) \n\\Dxu\\. 


2.2. Viscosity Solutions. Becanse we will be considering viscosity 
solutions, we will recall the main dehnitions and properties. We begin 
with the Riemannian jets and Jx~■ (See [Il[3] for a more complete 
analysis of such jets.) 

Definition 1. Let u be an upper semi-continuous function. Consider 
the set 

K^~^u{x) = |(p G (7^ in a neighborhood of x, p{x) = u{x), 

^{y) A u{y), y ^ X in a. neighborhood of xj-. 

Each function p G K‘^'^u{x) determines a vector-matrix pair ( 77 , X) via 
the relations 

T] = [Xiip{x),X2ip{x),... ,Xnp{x)) 

W, = i(W(V,(v.))(x)+V,(W((^))(x)). 
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We then define the second order superjet of u at x by 

= {{r],X) : ip e K‘^’'^u{x)}, 
the second order subjet of u at x by 

j‘^~u{x) = —j‘^’~^{—u){x) 
and the set-theoretic closure 

Jx^u{x) = {{r],X) : 3{xn,r]n, Xn}neN with G j‘^^u{xn) 

and {Xn,u{Xn),'nn,Xn) “t {x, u{x), 7], X)}. 

We then use these Riemannian jets to define viscosity cxD(a;)-harnionic 
functions as follows: 

Definition 2. A lower semi-continuous function v is viscosity cxd(x)- 

superharmonic in a bounded domain hi if n ^ cxd in each component 

_ 2 — 

of Q and for all Xq E fi, whenever (^,3^) G Jf n(xo), we have 

-((K.o + ll?f «.-DanpW>ln||e||) > 0. 

An upper semi-continuous function u is viscosity cxD(x)-subharnionic 
in a bounded domain hi if m ^ — cx) in each component of hi and for all 
Xq G Q, whenever {ri,X) G Jf u{xo), we have 

-((Ar/,r/) + ||r/||2(r/,D^lnp(a;))ln||r/||) < 0. 

A function is viscosity cxD(a:)-harmonic if it is both viscosity cxD(a:)- 
subharmonic and viscosity cxD(a;)-superharmonic. 

Similarly, we have the following definition concerning Ax,p(a;)'U. 

Definition 3. A lower semi-continuous function v is viscosity p(x)- 

superharmonic in a bounded domain hi if n ^ cxd in each component 

_ 2 — 

of Q and for all Xq G fl, whenever ((^, y) G Jf v{xo), we have 

An upper semi-continuous function u is viscosity p(a;)-subharnionic 
in a bounded domain hi if m ^ — cx) in each component of hi and for all 
Xo G fl, whenever {rj^X) G Jf m(xo), we have 

A function is viscosity p(x)-harmonic if it is both viscosity p(t)- 
subharmonic and viscosity p(x)-superharmonic. 
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Remark 1. In the above definitions, we may replace the right-hand 
side of each inequality by an arbitrary function. In that case, we use 
the term viscosity oo{x)-sub solution, etc. 

Our main tool is the Riemannian Maximum Principle [1], which we 
include for completeness. 


Theorem 2.1. Riemannian Maximum Principle Let u be upper 
semicontinuous in a bounded domain Vt C M"'. Let v be lower semicon- 
tinuous in O. Suppose that for x G dfl we have 


limsupM(|/) < hminfu(i/), 

y^x 

where both sides are not +cxo or —oo simultaneously. If u — v has a 
positive interior local maximum 


sup(u — u) > 0 
o 

then we have: 

For T > 0 we can find points Xr,yr G M"" such that 

i) 

lim Tfi{xr,yr) = 0, 

r—^oo 

where 

fi{x,y) = \x-y\^, 

for a fixed a > 2. (That is, is a power of the Euclidean 
distance.) 

ii) There exists a point h G O such that Xr ^ x (and so does y^- by 
(i)) and 

sup(u — v) = u{x) — v{x) > 0, 
n 

hi) there exist symmetric matrices Xr^yr and vectors pf, pI so 
that 


b) 

c) 

d) 


and 


{pt,Xr) G f^'^u{xr), 


pf -pf = o(l) 


as T ^ oo. 


Xr<yr + 0 ( 1 ) 
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We will also use the following corollary. 

Corollary 2.2. |3l Corollary 2.1] Let u and v be as in Theorem \2.1{ 
and in addition, let one of u or v be locally Lipschitz. Let a = 2. 
Then, for the vectors rjf and ri~ and the matrices AV and W as in the 
theorem, we have 

iiw^ir-iiwir = 0(1) 

and (A.,77+,77+) - (Whr )W) = o(l). 

This corollary is a consequence of the facts that the Riemannian 
distance is comparable to the Euclidean distance and that the choice of 
penalty function y) is the square of the Euclidean distance. Note, 
however, that even though the vectors pf and rj~ are not necessarily 
equal, we can still produce key estimates. 

3. Existence-Uniqueness of oo(x)-harmonic functions 

Let U be a bounded domain in M" and / : dfl —>■ M be a (Riemann¬ 
ian) Lipschitz function. 

We will hrst establish the existence of cxD(a;)-harmonic functions using 
Jensen’s auxiliary equations [5]: 

min{||L)xn|p-£, -A;:£^oo(x)m} = 0 and max{e-\\Dxu\\'^,-Ax,oo(x)u} = 0 

for a real parameter e > 0 . The procedure for existence of viscosity 
solutions to these equations (and viscosity cx)(a:)-harmonic functions) 
is identical to Section 4] and [HI Section 2], up to the obvious mod- 
ihcations. For completeness, we state the steps as one theorem and 
omit the proofs. 

Theorem 3.1. [SI 12] We have the following results: 

(1) Let e G M. If Uk is a continuous potential-theoretic weak sub- 

(super-)solution with u G to: 

J -Ax,kp(x)Uk = in Q 

\ u = f on dVt 

then it is a viscosity sub-(super-)solution. 

(2) Letting A: —)■ oo, we have Uk —)■ Uoo uniformly (possibly up to a 
subseguence) in U with u^o G 1 U^’°°(U) fl ^(U). 

(3) The function u^o is a viscosity solution to 

mindiDoenoolP - -Ax,oo(a:)Moo} = 0 when £>0 
max{£ - \\Dxu^\\^,-Ax^oo{x)Uoo] = 0 when e < 0 
-Ax,oo(x)Moo = 0 when e = 0. 
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In light of [21 Lemma 5.6] and [HI Lemma 2.2], the Main Theorem 
follows from showing the nniqueness of viscosity solutions to the Jensen 
auxiliary equations. We will establish this result, and point out where 
we digress from the Euclidean proof. 

Theorem 3.2. Let v = Uoo be the viscosity solution from Theorem \S.l\ 
to 

min{||DxM|p - e, -A^,oo(x)R} = 0 (3.1) 

in a bounded domain If u is an upper semi-continuous viscosity 
subsolution to Equation fl3.ll) in so that u < v on dfl, then u < v in 

n. 

Proof. Following [HI Lemma 3.1] and [21 Theorem 5.3], we may assume 
WLOG that n is a strict viscosity supersolution. Suppose 

sup(ti — n) > 0 
o 

and let '0(x,j/) = |x — y\‘^ be the square of the Euclidean distance 
between the points x and y. By the Riemannian Maximum Principle 
(Theorem 12.ID . there are points Xr and Pr in (for sufficiently large r) 
with the property that there are vectors pf., pf and symmetric matrices 
Xr , y-r SO that 

{pf, Xr) e u{Xr) and {pf,yr)e.f’ v{yr). 

Since n is a viscosity subsolution and v a strict viscosity supersolu¬ 
tion, we have, for some p > 0, 

0 > mm{\\pff - e,-{Xrpf,pf)-\\pf\\^{pf,Dxlnp{xr))ln\\pf\\} 

0<p < mm{\\pff-e,-{yrP;,pf)-\\p;f{p;,Dx\np{yr))\n\\p;\\}. 

Subtracting these equations, we obtain 

0 <^ < max{\\pf\\‘^-\\pff,{Xrpf,pf) - {yrPf,p;) 

+ \\vtf{vt,Dxlnp{xr))ln\\pf\\ (3.2) 

- Ihr f (W , Dxlnp{yr)) In \\p;\\}. 

Here is where the proof diverges from the Euclidean case. In the 
Euclidean case, the vectors pf and pf are equal, rapidly leading to a 
contradiction in Equation fl3.2p as r —)■ cxd. However, in the Riemannian 
environment, these vectors are not, in general, equal. So, we will have 
to estimate the right-hand side more carefully. 

Since v is locally Lipschitz, can invoke Corollary 12.21 to obtain 

WVyrf-WVxrf^O and {Xrpf ,pf) - {yrPf ,pf) ^ 0 


as r —)■ cx). 
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Thus, we are left to consider 

Dxlnp(x^)) In \\r]^\\ - \\r];f{r];,Dx Inp(y^)) In \\t];\\. 
We begin by expressing the sum as 

Wvtfivt, Dxln^ixr)) In \\7]+\\ - \\r]-f{r];, Dx lnp( 2 /^)) In \\r]; 
= \\vtf{vt,Dxh-i^{xr))(^ln\\r]^\\ -ln||77-||^ 

+ (^llhr IP - IIW lP^(^P,^xlnp(a:r))ln|h-|| 

+ WVrf ({vt^ Dxlnp{xr)) - {r];, Dxln^ixr))^ ln||?7"|| 

+ \\Vrf({Vr^Dx\np{xr)) - (t]", Inp(2/r))^ ln||w II 

= I + II + III + IV. 

To estimate each term, we will need the following Lemma. 


Lemma 3.3. For some constant K > 0, the vector 7].^. satisfies 

\fe < ||?7“|| < K. 

In addition, for that same K, and sufficiently large r, the vector rj'^ 
satisfies 

^^/i< ||r7+|| < K. 

As a consequence, there is a constant L so that 


\n\\r]; 


< L < oo. 


Proof. Since n is a strict supersolution to Equation fld.lh . we have 0 < 
e < II^^Llp. Next, since v is locally Lipschitz, the proof of Corollary 12.21 
shows there is a constant C so that 

Tf){Xr,yr)^ < C. 

By the proof of the Riemannian Maximum Principle fTheorem 12. ip . 
gf = -TA{y.r)Dy{fi{Xr,yr)) 

where Dy denotes Euclidean differentiation in the y-variable and A(?/t-) 
is the coefficient matrix of the frame at y^ in terms of the standard 
Euclidean vectors (Section 2.1). Because fi^Xr^yr) = \xr — 2/rP and 
yr & It, Si bounded domain, we conclude that for some hnite constant 
C'l, 
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Similarly, 


7 ^+ = Tk{Xr)D^{'ll){Xr,yr)) 


SO that for some finite constant 6 * 2 , 

llh+ll <C'2. 

By part iii(c) of the Riemannian Maximum Principle (Theorem I2.ip . 
for sufficiently large r, 0 < \\/s < ||77^||. The Lemma then follows. □ 

Term I: The absolute value of the hrst term is controlled by 


WVrfWDx lnp(a;^)||Lo 


In II77+II - In II 77 ; 


Using Lemma 13731 and the fact that 1 < p(a:) G fl we 

have that this is, in turn, controlled by 

C'(^ln||77+|| -ln||77; 

for some hnite constant C. Suppose that 

Wvt 


Then 


ln||77+|| -ln||77; 


In II 77 +II - In II 77 ; 


In 


\\V~r 


< In ( 1 + 


= In 


Wvt -V~r 




\\V~r 


\\v; 


Lemma 13.31 and the Riemannian Maximum Principle then imply as 
r —)■ cxo, 

^ 0 . 


If instead. 

In II 77 +II - In II 77 ; 


In II 77 +II - In II 77 , 

II 


In 


11 ^; 


= — In 


11 ^^ 


II??; 


= In 


WV'r 


WV^r 


0 


then a symmetric argument gives 

In II 77 +II - In II 77 ; 

as r —)■ 00 , so that 

llhr ir(??r ,^Tlnp(x^))(^ln||77+||-ln||77-|| ) ^0 

as r —)■ 00 . 

Term II: By Corollary 12.21 
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Using Lemma ESI we have 




< C'||Dxlnp(x^)||Loo. 


As in Term I, this term is bounded and so as r ^ cxo 

- \\Vrf'^{vt:Dxlnp{xr))ln\\r];\\ 0 . 

Term III: As in the previous terms, the absolute value of this term is 
controlled by 

C\\vt -Vrl 

By the Riemannian Maximum Principle, we then have as r —)■ oo, 

\\VTf({vt,Dx\np{xr)) - {rj;, Dx\np{xr))^ ln||?7"|| ^ 0. 

Term IV: By the Riemannian Maximum Principle, both Xr and y-r 
converge to a point x. By the regularity of p(x), we have the absolute 
value of Term IV is controlled by 

C'fT':rlnp(x^) - L>;:£lnp(?/^)^ -^0 


as r —)■ cxD. 


□ 


An analogous argument produces the following Corollary. 

Corollary 3.4. Let v = Uoo be the viscosity solution from Theorem \3.1\ 
to 

max{e - ||T>xm|P, -Ax,oo(a;)M} = 0 (3.3) 

in a bounded domain fi. If u is an lower semi-continuous viscosity 
supersolution to Equation fl3.3p in so that u > v on dQ, then u > v 
in f2. 


4. A Harnack Inequality 

We include a Harnack inequality for completeness. First, we have 
the following lemma whose proof is identical to [HI Lemma 4.1] and 
omitted. 


Lemma 4.1. Let u be a positive viscosity oo{x)-harmonic function and 
C a positive, compactly supported smooth function. Then 


sup 


DxC{x)Dx \nu{x) 


p{x) 

< sup 


DxC{x)+C{x) In 



Dx Inp(x) 


p{x) 


As in [6l Section 4], we have the following Harnack inequality as a 
consequence. 
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Theorem 4.2. Let u be a positive viscosity oo{x)-harmonic function. 
Let Bj. be a ball of radius r > 0 contained in the bounded domain 
Let B 2 r be the concentric ball of twice the radius also contained in 
Then 



for some constant C depending on u{x). 
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